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ABSTRACT 

Giant H 11 regions (GHIIRs) in nearby galaxies are a local sample in which we can study in 
detail processes in the interaction of gas, dust, and newly formed stars which are analagous 
to those which occurred in episodes of higher intensity in which much of the current stel- 
lar population was born. Here, we present an analysis of NGC 588, a GHIIR in M33, 
based on optical Integral Field Spectroscopy (IFS) data obtained with the PMAS instru- 
ment at the 3.5 m telescope of Calar Alto Observatory, CAHA, together with Spitzer in- 
frared images at 8 /im and 24 /im. The extinction distribution measured in the optical shows 
complex structure, with three maxima which correlate in position with those of the emis- 
sion at 24 /im and 8 fim. Furthermore, the Ha luminosity absorbed by the dust within 
the H II region reproduces the structure observed in the 24 image, supporting the use 
of the 24 /im band as a valid tracer of recent star formation. A velocity difference of 
^50 km s^^ was measured between the areas of high and low surface brightness, which 
would be expected if NGC 588 were an evolved GHIIR. We have carefully identified the 
areas which contribute most to the line ratios measured in the integr ated spe c trum. Those 
line ratios which are used in diagnostic diagrams proposed by Baldw iiTet al.l (Il98l[) show 
a larger range of variation in the low surface brightness areas. The ranges are ^0.5 to 
1.2 dex for [N ll] A6584/Ha, 0.7 to 1.7 dex for [S ll]AA6717,6731/Ha, and 0.3 to 0.5 dex 
for [O III] A5007/H;3, with higher values of [N ll] A6584/Ha and [S ll] AA6717,6731/Ha, and 
lower values of [O ill] A5007/H/3 in the areas of lower surface brightness. Ratios correspond- 
ing to large ionization parameter ([/) are found between the peak of the emission in H/3 
and the main ionizing source decreasing radially outwards within the region. Differences be- 
tween the integrated and local values of the U tracers can be as high as ^0.8 dex, notably 
when using [O ill] AA4959,5007/[O ll] AA3726,3729 and in the high surface brightness spax- 
els. [Oll]AA3726,3729/H/3 and [Olll]AA4959,5007/[Oll]AA3726,3729 yield similar local 
values for the ionization parameter, which are consistent with those expected from the inte- 
grated spectrum of an H II region ionized by a single star. The ratio [S ll] AA6717,6731/Ha 
departs significantly from the range predicted by this scenario, indicating the complex ioniza- 
tion structure in GHIIRs. There is a significant scatter in derivations of the metallicity using 
strong line tracers as a function of position, caused by variations in the degree of ionization. 
The scatter is smaller for N203 which points to this tracer as a better metallicity tracer than 
N2. One interesting result emerges from our comparison between integrated and local line 
ratio values: measurements of the line ratios of GHIIR in galaxies at distances >25 Mpc may 
be dominated by the ionization conditions in their low surface brightness areas. 

Key words: H ll regions: individual: NGC 588 - galaxies: individual: M33 - stars: Wolf- 
Rayet - ISM: abundances - ISM: kinematics and dynamics - dust, extinction. 
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1 INTRODUCTION 

Large areas of ionized gas known as Giant H II regions (GHI- 
IRs) constitute the mos t conspicuous places of star formation 
in normal galaxies (see IShieldsl [l990l. for a review). Their di- 
ameters typically range between ~100 p c to ~800 pc (e.g. 
lKennicutliri984lAlonso-Herrero et al.ll2002l) while their Ha lumi- 
nosity range expa nds up to 3 orders of magnitudes (^ 10'^^ — 
10"° erg s~\ e.g. iKennicutll [19 84: Rozas et al. 199l| iFirpo et al.l 
|2005| ; iMonreal-Ibero et al J 120071 : IRelano & Kennicutj 120091) . Re- 
garding their morphologies, some of them present a compact dis- 
tribution with high surface brightness while others have a more dif- 
fuse emission. Also, they can present multiple cores and/or shells 
or ring-like features. In the same manner, they present very dif- 
ferent content of gas, varying between ^ 10'' Mq to almost 
10^ Mq while they usu ally have ~ 10^ — 10^ M© in stars (e.g. 
ICastellanos et al.ll2002bl) . Finally, GHIIRs are small-scale exam- 
ples of the extreme ev ents of star formation occurring in star- 
burst galaxies (e.g. Alon so-Herrero e t al. 2009; Garcfa-M arm et all 
l2009l) . Thus, a good knowledge of these objects is highly valuable 
for a better understanding of these more violent phenomena. 

The variety of properties in GHIIRs implies a complex struc- 
ture, far from the textbook-like Stromgren sphere. Moreover, mod- 
elling of star-forming regions showed that in case of multiple 
ionizing sources, geometrical effects affect the physical proper- 
ties (i. e electron tem perature and ionization structure) of these re- 
gions tercolano et a l. 2007; Jamet & Morisset 2008). The inhomo- 
geneities of the Interstellar Medium (ISM) have been taken into 
account within the Hll regions using the filling factor. This de- 
scribes the fraction of the total volume of the H II region with high 
dense gas while the remaining volume is considered to be of negli- 
gible density ( Osterbrock & Flather 1959). Recently, detailed mod- 
els show that the density variations assuming optically thick high 
density gas clumps give ri se to inhomogeneities in the tem perature 
and ionization parameter l lGiammanco et al.ll2004l l2005h . Thus, a 
single value per physical magnitude, usually extracted from a spe- 
cific area of the region, is not necessarily representative of its phys- 
ical conditions. 

It is in very nearby GHIIRs (i.e. at Z)^6 Mpc), where ground 
based optical telescopes can achieve high linear spatial resolution 
(i.e. ^40 pc arcsec"^), enough to resolve the different elements 
(i.e. star clusters, ionized gas, dust, etc.), playing a role in the inter- 
action between the massive stars and its surrounding environment. 
Also, it is there where the variations of the physical and chemical 
properties of the ionized gas can be properly sampled. 

From the observational point of view, the study of these re- 
gions would benefit from high quality spectroscopic data that map 
in an un-biased way the surface of the GHIIR. Nowadays, the 
technique of integral field spectroscopy (IFS), able to record si- 
multaneously the spectra of an extended continuous field, offers 
the possibility of performing such a mapping. At present, studies 
of GHIIRs based on IFS are still scarce. An example is provided 
by [Gar cfa-Benito et al. (2010), who anal yzed one of thebright- 
est GHIIRs in NGC 6946. Also. iLopez-'sanchez et al.l ( l2010h pre- 
sented a detailed study of a star-forming region at the lower limit 
of GHIIRs in terms of Ha luminosity and size in IC 10, our clos- 
est starburst. Recently, we presented a 2D spectroscopic analysis of 
the second brightest H II regio n in M 33, the Triangulum Galaxy: 
NGC 595 ( lRelanoetal]|2010h . There, we showed how the optical 
extinction map and the absorbed Ha luminosity are spatially cor- 
related with the 24^m emission from Spitzer, and how the ioniza- 
tion structure of the region nicely follows the Ha shell morphol- 



ogy. Moreover, we evaluated the reliability of different line ratios 
as metallicity (Z) tracers. In a companion paper, we presented a 
novel approach to model these complex structures. There, we re- 
produced our observations by jointly fitting the radial profiles of 
different optical (i.e. Ha, [O II] AA3726,3729, [OlIl]A5007) and 
infrared (i.e. 8 nm, 24 ^i m) magnitudes to a set of C LOUDY-based 
photoionization models jPerez-Montero et"al]|20IOl) . 

Our experience with NGC 595 shows the importance of carry- 
ing out an analysis using the combined information of optical and 
infrared data together with modelling. However, given their diver- 
sity, the sample of GHIIRs studied by means of this methodology 
cannot be reduced to only one example. Instead, studies of other 
regions sampling a different range in the parameter space would be 
desirable. We present here, the analysis of a region with a different 
morphology, and relatively lower metallicity and high Ha luminos- 
ity: NGC 588. This region is located in the outskirts of M 33, at the 
end of a spiral arm, at a radius of ~I4' (i.e. ~ 3.42 kpc for a dis- 
tance to M33 of 840 kpc;|F reedman et aljl99lh. Thi s area presents 
H I emission dVerlev et all201Q ; ,Gratier et al.lf20IOh . However, no 



local H2 (or CO) emission has been detected towards NGC 588 
JVerlev et al.ll2010l ; llsrael et aljl990h . With a size of 30" x 50" 
(i.e.~120 pcx200 pc at our assume d distance), NGC 588 has been 
classified within the ring-like class jSabbadin et alj|l980l) . Its stel- 
lar content has been thoroughl y studied 1 Jamet et alj2004l ; [Pellerinl 
l2006l ; IUbeda & Drisserj|2009h . The different estimates for its total 
stellar mass range between ~ 1.3 x 10'' and ~ 5.6 x 10'' M0 with 
an age for the burst of ~ 3.5 — 4.2 Myr and low metallicity (i.e. 
Z ~ OAZq), consist ent with direct meas urements using long-slit 
(12+log(0/H) = 8.30, lvnchez et alj|l988l) . NGC 588 contains two 
Wolf-Rayet (WR) stars. One of them was classified as WNL while 
the other as Of peAVN9, an intermediate object between Of and 
WN stars ( Massev et al ][l996i). From the kinematic point of view, 
TAURUS-2 Fabry-Perot data show how NGC 588 seems a rela- 
tively evolved system, with its ionized g as kinematics dominated b y 
a collection of large stellar wind shells jMufioz-Tufion et al .119961) . 
Finally, a relatively faint and point-like X- ray emitting source as - 
sociated with this region has been detected jPlucinskv et al.ll2008l) . 

In this work, we combine 2D optical spectroscopic observa- 
tions with the Potsdam Multi-Aperture Spectrophotometer (PMAS, 
fRoth et al. 2005) and infrared imaging with Spitzer. These will con- 
stitute the first published results from data obtained with the new 
PMAS's CCD. A modelling of NGC 588 intending to reproduce 
the observed magnitudes will be presented in a companion paper 
(Perez-Montero et al. in preparation). We describe the observations 
and data processing to create the maps of the relevant magnitudes in 
Section |2] The main observational results are presented in Section 
[3] Finally, we summarize our main conclusions in Section|4] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The IFS data of NGC 588 were obtained on October 9-10, 2009 
during the commissioning r un of the new Pots dam Multi-Aperture 
Spectrophotometer (PMAS. iRoth et aljlioosh CCD at the 3.5 m 
telescope at the Calar Alto Observatory (Spain). The new PMAS 4k 
x 4k CCD is read out in fo ur quadrants which have slightly differ- 
ent gains ( Roth et alfcOIOl) . Data were taken using the Lens Array 
Mode (LARR) configuration which is made out of a 16 x 16 array of 
microlenses coupled with fibres (hereafter spaxels). We used the 1" 
magnification which provides a field of view (FoV) of 16" x 16". 
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Figure 1. Mosaic to map NGC 588 overplotted on a continuum-subtr acted 
Ho direct image from NOAO Science Archive jMassev et aDl2007h . Tlie 
orientation is north up and east to the left. The Ha image is shown in log- 
arithmic stretch to better enhance all the morphological features of the H II 
region and covers a range of 3.4 dex. 



We used the V600 grating and the 2x2 binning mode achieving 
an effective dispersion of 1.59 A pix~^, and a ~3.4 A full width 
half maximum spectral resolution. With the new PMAS's CCD, the 
covered spectral range was from 3 620 to 6 800 A for most of the 
spaxels. This permits us to observe the main emission lines in the 
optical from [O II]AA3726,3729 to [S II]AA6717,6731. However, 
for ~40 spaxels, always at the edge of the LARR, this spectral 
range was slightly reduced due to the vignetting associated with 
the 3.5 m telescope. This prevents us from obtaining information 
for the [Oll]AA3726,3729, [S II]AA6717,6731, and Ha emission 
lines in some specific areas (see also section|2j2}. 

We made a mosaic of 6 tiles to map most of the surface of 
NGC 588. The distribution of the different tiles is shown in Fig.[T] 
overplotted on the Ha emission-line image from Na tional Optical 
Astr onomy Observatory (NOAO) Science Archive jMassev et al.l 
l2007i) . Contiguous tiles had a 2'.'0 overlapping, to make easier a 
common relative flux calibration of the data. In total, we covered a 
field of 30'.' Ox 44'.' which at the distance of NGC 588 corresponds 
to ~ 120 pcxl80 pc. 

We obtained three exposures of 400 s per tile. Atmospheric 
conditions during the observations were non-photometric and typi- 
cal seeing ranged between l'.'2 and l'.'6. The sky transparency dur- 
ing the observing nights presented variations of ^15 %. All the data 
were taken at airmasses ^1.1 in order to prevent strong effects due 
to differential atmospheric refraction. 

In addition to the science frames, continuum and HgNe arc 
lamp exposures in order to minimize the effects due to instrument 
flexures. Also, a nearby sky background frame was obtained during 
the second night by moving the IFU off-target. 

Finally, exposures of the spectrophotometric standard star 
BD-I-28D4211 were obtained in order to correct for the instrument 
response and perform a relative flux calibration. 



2.2 Data reduction and map creation 

The first steps of the data reduction were done through the P3d 
tool that is designed to be u sed with fibre-fed integral-field spectro- 
graphs jSandin et al.ll2O10h . After trimmimg, combining the four 
quadrants and subtracting the bias level, the expected locations of 



the spectra were traced on a continuum-lamp exposure obtained 
before each target exposure. We extracted the target spectra by 
adding the signal from the central traced pixels ant its four neigh- 
bours. The spectra were wavelength calibrated using the exposures 
of HgNe arc lamps obtained immediately after the science expo- 
sures. We checked the accuracy of the wavelength calibration us- 
ing the [Ol]A5577 A sky line, and found standard deviations of 
<0.14A, which allowed us to determine the centroid of line with 
an accuracy of ~7 km s^^ at ~5 000 A. 

Fibres have different transmissions that may depend on the 
wavelength. The continuum-lamp exposures were used to deter- 
mine the differences in the fibre-to-fibre transmission and to ob- 
tain a normalized fibre-flat image, including the wavelength depen- 
dence. This step was carried o ut by running t he f iber.f lat . pi 
script from the R3D package jSanchezll2006h . In order to homoge- 
nize the response of all the fibres, we divided our wavelength cal- 
ibrated science images by the normalized fibre-flat. To estimate 
the accuracy of the fibre-to-fibre response correction, we fitted a 
Gaussian to four emission lines distributed along the whole spectral 
range in an extracted, wavelength calibrated and flat-field corrected 
arc exposure. We used the ratio between the standard deviation and 
the mean flux in each line as a proxy for the accuracy of the fibre-to- 
fibre response correction. For those lines in the central part of our 
spectral range, this correction was very good, with ratios of ~2%. 
In the blue and red edges these ratios reached values of ~12%, due 
to the contribution of fibres affected by vignetting. 

In the next step, the three exposures taken for the same point- 
ing were combined in order to remove cosmic rays, using the 
imcombine routine in IRAfQfIux calibration was performed us- 
ing the IRAF tasks standard, sens fun c and calibrate. We 
co-added the spectra of the central fibres of the standard star expo- 
sure to create a one dimensional spectrum that was used to obtain 
the sensitivity function. For wavelengths larger than 4 000 A, the 
uncertainty in the flux calibration is ~ 1%, while across the bluer 
spectral range (i.e. <4 000 A), the associated error can reach ~ 5%. 

Given the large size of the new PMAS's CCD, th e four cor- 
ners o f each exposure suffer from telescope vignetting jRoth et al] 
I2OIOI) . Due to the way the fibres of the LARR are arranged at the 
entrance of the spectrograph, a maximum of two columns of fibres 
at the east and west sides of the LARR were affected by this. The 
wavelength range affected was larger for those fibres located more 
towards the edge of the spectrograph and could be up to ~3 898 A 
in the blue end and from ~6 498 A in the red end. These parts of 
the spectra were masked and then we used the offsets commanded 
to the telescope and the PMAS acquisition images to construct a 
mosaic datacube. 

After creating this datacube, maps for the different ob- 
servables we r e deri ved following the methodology presented in 
iRelaiio et alj ( I2OI0I) . Basically, we performed a Gaussian fit to 
the emission line s using the IDL-based routine mpfitexpr 
jMarkward^ l2009h and derived the quantities of interest for each 
individual spaxel. Then, we used these together with the position 
of the spaxels within the datacube to create an image suitable to 
be manipulated with standard astronomical software. Hereafter, we 
will refer to this with both terms: map and image. 

Even after having allowed for 2" of overlap between tiles, vi- 



^ IRAF is distributed by the National Optical Astronomical Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Table 1. Observed and extinction-corrected emission line fluxes, normal- 
ized to I(H/3)=1 from the integrated spe ctrum of NGC 588 . f(A) is the red- 
dening function normalised to H/3 from lFluks et al.l jl994l) . 



Emission Line 




f(A) 


I(A)/l(H/3) 


I(A)/I(H/3)corr 


3727 [Oil] 


3725.52 


0.255 


1.922±0.048 


2.317 ± 0.091 


3868 [Nclll] 


3866.01 


0.227 


0.389±0.013 


0.461 ± 0.020 


3889 HS+Hcl 


3886.21 


0.223 


0.423±0.014 


0.501 ± 0.021 


3970 He+[Ne III] 


3966.03 


0.208 


0.152±0.011 


0.179 ± 0.013 


4101 Hc5 


4099.02 


0.180 


0.199±0.0I3 


0.274 ± 0.016 


4340 H7 


4338.19 


0.133 


0.397±0.010 


0.451 ± 0.020 


4861 


4859.46 


0.000 


l.000±0.016 


1.000 ± 0.021 


4959 [0 III] 


4957.02 


-0.027 


l.500±0.021 


1.524 ± 0.046 


5007 [Olll] 


5004.90 


-0.041 


4.400±0.049 


4.428 ±0.131 


5876 He I 


5872.40 


-0.217 


0.154±0.004 


0.140 ± 0.007 


6563 Ha 


6559.34 


-0.314 


3.508±0.034 


2.938 ± 0.114 


6584 [Nil] 


6579.40 


-0.316 


0.246±0.004 


0.210 ± 0.008 


6678 He I 


6674.29 


-0.329 


0.036±0.002 


0.031 ± 0.002 


6717 [S II] 


6712.40 


-0.332 


0.210±0.004 


0.178 ± 0.007 


6731 [Sll] 


6726.76 


-0.335 


0.147±0.003 


0.124 ± 0.005 



Table 2. Main diagnostic emission line ratios and integrated physical prop- 
erties of NGC 588. All the physical properties were estimated using our 
extinction corrected integrated fluxes. 



Parameter 


ValuecoiT 


log([NlI]A6584/Ha) 


-1.146±0.077 


log ([S II]AA6717,6731/Ha) 


-0.988±0.078 


log([Olll]A5007/H/3) 


0.646±0.05() 


log ([OIII]AA4959.50()7/[OII]A3727) 


0.406±0.069 


log R23 


0.915±0.053 


N203 


-1.792±0.127 


[Sll]A6717/[Sll]A6731 


1.435±0.143 


log([N II] A6584/[0 II] AA3726,3729) 


-1.043±0.077 


log([N II] A6584/[S 11] AA6717,6731) 


-0.271 ±0.078 


c(H/3) 


0.26±0.03 


Ue (cm~3) 


< 100 


qeff (cm s-1) 


3.9 X lO'' 


Te (K) 


11 140±180(*) 


12-l-log(0/H) 


8.16±0.02 



WValue derived from the [OlII] A4363/A5007 line ratio bv I Jamet et al] 
i2005l) . 

gnetting prevents us from deriving information for some lines in 
the central columns of our mosaic. In those maps where (at least) 
one line affected by vignetting was involved, those spaxels were 
masked and interpolated using f ixpix within IRAF, for presenta- 
tion purposes. 



3 RESULTS 

3.1 Integrated properties 

Here we analyze the integrated spectrum for NGC 588, obtained 
after co-adding the signal of ~ 1 000 spaxels. The sky substracted 
spectrum is shown in Fig.|2l where we used two different normal- 
ization factors in order to better display all the observed emis- 
sion lines. The positions of the detected nebular emission lines 
were marked with labels. In addition, several sky line residuals 
are clearly visible. Flux for the main emission lines was mea- 
sured using splot within IRAF, which integrates the line inten- 



sity ov er a locally fitted continuum. As was shown in lRelafio et al.l 
l l2010l) . both splot and mpf itexpr give similar results for the 
high signal-to-noise integrated spectrum. We derived the redden- 
ing coefficient, c(H/3), from the F(H(5)/F(H/3), F(H7)/F(H/3) and 
F(Ha)/F(H/3) line ratios. We performed a linear fit (by minimiz- 
ing the chi-square error statistic) to the difference between the the- 
oretical and obser ved B aimer decrements vs. the reddening law 
l lFIuksetaI.lll994h . while simultaneously solving for the effects 
of underlying Balmer absorption with equivalent width, EWats. 
We assumed that the EWg ba is the same for all Balmer lines 
(e.g. i Kobulnickv et"Zlll999l). The theoretical Balmer 1 ne intensi- 



Storev & Hummeil i ll 9951) assuming Case 



ties were obtained from 

B,Te = 10"^ K, ne=100 cm (typical values found in H II regions, 
lOsterbrock & Ferlandll2006h . 

We derived a mean extinction for the region of Av = 0.548. 
Given that the Galactic extinction in the line of sight to NGC 588 
is Ay"^^ = 0.146 (Schlegel et al. 1998) most of it is intrinsic to 
NGC 588 and agrees within the uncertainties with previous mea- 
surement on its brightest parts {Ay = (1 49 - 0.81. lMelni"c3ll979l: 
IViallefond & Goss|]l986l : iMelnick et al ]fl987l) . 

Table [T] presents both the measured and extinction-corrected 
fluxes. Errors were est imated using the formula presented in 
ICastellanos et al.l jiooial) . 



1 /-) EW 

O-lino = 0-cont X iV'/'(l + -^^^ ) 



,1/2 



(1) 



where tJcont is the standard deviation in a continuum close to the 
line of interest. A*' is the number of pixels sampling the line, EW 
is its equivalent width and A A is the dispersion in A pix~i. The 
quoted uncertainties of the reddening corrected line fluxes take into 
account the measurement and reddening errors. 

To investigate if the nebular properties derived from long-slit 
are representative of the whole H II region, we compared our flux 
measurements from the integrated spectru m of NGC 588 to the 
values presented by Jamet et al. (2005) and IVflchez et al.) h988l) 
who used long slit at position angle (RA.) ~ —45° i ncluding the 
main ionizing cluster in NGC 588 (e.g. see Fig. I in I Jamet et ahl 
I2OO5I) . Most of our measurements present a difference with respect 
to H/3 when comparing with previously reported measurements that 
can range between ~8% (e.g. [OlIl]A5007 and HelA6678) up to 
~64% in the case of H8-l-He IA3889. This indicates that it is not 
trivial how the line ratios measured at the brightest knots trace those 
for the H II region as a whole. This issue will be explored in more 
detail in Sec. 13. 61 

Table [2] contains the principal diagnostic emission-line ratios 
measured from the integrated spectrum as well as the physical pa- 
rameters derived from them. Given the blue-shift for IVI 33 and the 
presence in the sky substracted spectrum of strong residuals for the 
HglA4358 sky-line, it was not possible to measure the [O III] A4363 
nebular line and thus determine the electron temperature (Te) by 
means of the [O III] A4363/A5007 line ratio. Another possibility 
would have been using [Sll]A4067 and/or [Nll]A5755 together 
with [Sll]AA6717,6731 and/or [Nll]A6584. To search for these 
lines, in addition to the total spectrum, we created a spectrum by 
co-adding the spectra in an aperture of 5 x 9 spaxels centred at the 
peak of emission for the ionized gas. In this way, we increased the 
signal-to-noise ratio - and thus, our detection limit - since only 
the brightest spectra were included. However, these features were 
detected neither in the total spectrum nor in the one involving the 
brightest spaxels. T hus, we assumed a of 1 1 140 K, as derived by 
Ijamet et alj j2005l) from the [O III] A4363/A5007 line ratio to esti- 
mate the electron density (jie) and used the task temden, based on 
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Figure 2. Integrated spectram of NGC 588 obtained by co-adding the signal of all the spaxels in the field of view with two different normalizations to better 
visualize the different observed emission lines. Fluxes are in arbitrary units. The sky spectrum was created by combining the spectra corresponding to the 
vigneting-free fibres of the background frame. 



the f ivel program l lShaw&Dufouilll995 h included in the IRAF 
package nebular. The derived rie was consistent with being be- 
low the low density limit. 

Another quantity quoted in Table|2]is the ionization parameter, 
defined as: 



QjHo) 



(2) 



where Q{Ho) is the number of ionizing photons per second emitted 
by the stars, Rs is the Stromgren radius of the H II region and is 
the electron density. We estimated the number of ionizing photons 
using the expression provided by Kennicutt ( 1998,): 



Q{Ho){s'^) = 7.31 X 10" L(iyQ) (ergs" 



(3) 



and the extinction corr ected Ha luminosity reported by 
iRelano & Kermicut ^( 12009' ). Since the sulfur line ratio is consistent 
with being below the low density limit regime, we assumed a face 
value of Tie — 20 cm~'^. Also, we utilized a Stromgren radius of 
80 pc which is an approximated value inferred from the Hq image. 

Table |2] also includes several metallicity sensitive line ratios. 
The most widely used is probably the_R23=( [On] A A3726,3729-l- 
[Oiii]AA4959,5007)/H/3 index jPagel et alJ[T97^ . However, the 
Z — R2S relation is two valued and thus, independent metallic- 
ity tracers are needed to determine which of the two branches 
is appropiate for NGC 588. One possibility would be the 

= log([Nll]A6584/Ha) li ne ratio. According to the empiri - 
cal parametrization proposed bv IPerez-Montero & ContiniH2009l) . 
we derived a metallicity of 12 + \og(0/H) = 8.16. This re- 
sult does not point in a conclusive manner to either the up- 



per or the lower branch. Instead, the metallicity of NGC 588 
falls in the knee region of the Z ~ i?23 relation, where un- 
certainties can be as high as 0.7 dex. Alternatively, one can 
use the N203=log(([N II] A6584/Ha)/([0 III] A5007/H/3)). We de- 
rived a metallicity of 12 + logfO/ H) = 8.15 and 8.18 us- 
ing the parametrization p roposed by IPettini & Pagel ( l2004l ) and 
IPerez-Montero & Continil (j2009), respectively. Thus for the pur- 
pose of this section, we will consider as the characteristic metal- 
licity of the region, the mean of those derived from the N2 and 
N203 parameters: 8.16±0.02. This value agree within the uncer- 
tainties with the 12 + \og{0/H) = 8.17 metallicity reported by 
Jamet et al. (2005), is slightly lower than the expected value as- 
suming that this region follows the metall icity grandient for M 33 
(8.28±0.08, iRosolowskv & SimonI l2008h and ~0.15 dex lower 
than the value reported bv lvflchez et al.l l ll988l) . 



3.2 Structure of the ionized gas and the stellar component 

In Figure [3] we present the H/3 and [Oin]A5007 flux maps for 
NGC 588 derived from our PMAS data. Contours reproducing 
archive continuum images in the red and blue spectral bands have 
been overplotted for reference (see caption of Fig. [3] for details). 
Our PMAS data cover the whole southern part of the region plus 
most of the northern one. This GHIIR is dominated by emission 
from a broken elongated ring-like structure with major and minor 
axes of ~ 40" and ^ 25" (i.e. ^^160 pc and ~100 pc), respectively 
and at P.A.~10° and ~100°. In addition, there is a bridge of ion- 
ized gas emission joining the two halves of the ring-like structure 
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Figure 3. Left: Map of the observed H/3 flux derived from our PMAS data. Each spaxel has a 1" X 1" size. A logarithmic stretch covering 2.4 d ex was used to 
better enhance all the morphological features. Units are arbitrary. Contours correspond to the R band direct image from NOAO Science Archive jMassev et alj 
l2007h and show the location of the ionizing stars within the region. The orientation is north up and east to the left. The peak in this continuum image, at 
coordinates RA(J2000): Ih32m45.7s, Dec.(J2000): +30d38m55.1s, marks the origin of our coordinate system and will appear as a cross in the following 
figures for reference. Right: Similar map for the observed [O III] A5007 flux. The logarithmic stretch covers 3.3 dex and contours represent the HST-WFPC2 
image with the filter (program 5384; P.I. W. William), convolved with a (X'S Gaussian filter. 



from ~[1'.'0,~7'.'0] to ~[-7'.'0,-3'.'of| at P.A.~ -70°. The mor- 
phologies of the H/3 and the [OlIl]A5007 maps are very similar 
but show some rather subtle differences which indicate the com- 
plex ionization structure, which will be explored in more detail in 
section [X6l The main ionizing cluster, as depicted by the archive 
continuum images, is not at the centre of the ring like structure, but 
at ~2'.'0 from the peak of emission in H/3. Finally, there are also 
several secondary peaks of emission, most of them in the north- 
em half of the region, which are associated with very massive (i.e. 
30 — 45 M0) individual ionizing stars dJamet et al.ii2004) . 



3.3 Extinction distribution and dust 

The distribution of the extinction was derived by means of the Ha 
and H/3 emission line maps. We assumed an intrinsic Balmer emis- 
sion line ratio of Ha/H/3= 2.86 jOsterbrock & Ferlandl2006h for a 
case B aproximation and T e = 11 150 K and used the extinction 
curve of iFluks et al. I(ll994!). We included a 1 A correction to take 
into account the Yip absorption line due the underlying stellar pop- 
ulation. This absorption feature was clearly visible in a co-added 
spectrum extracted in an rectangular area with low surface bright- 
ness in the emission lines located at ~ [— 12'.'0, — 17'.'0] of about 



^ Hereafter, the reported positions will refer the relative coordinates to the 
main ionizing cluster. 



10" X 8". However, it was not detected, and thus impossible to be 
fitted for individual spaxels. 

The reddening map was cr eated assuming E{B — V) = 
Av/i.l teieke&Lebofskvl 19851) and is displayed in the left panel 
of Fig. |4] with the 24 ^m image from Spitzer overplotted with con- 
tours. This map shows how irregular the extinction distribution is 
with low values of redden ing (~ 0.00 — 0.2 5). This strongly con- 
trasts with the findings of Ijamet et al.l ( |2004|) who, using long-slit, 
reported an almost constant extinction of — — 0.11±0.02 
and reinforces the need of 2D unbiased spectral mapping to char- 
acterize the physical properties of GHIIRs. The optical reddening 
map presents three maxima that spatially correlate very well with 
the maxima of dust emission in the Spitzer 24 jim and 8 fim bands 
(see Fig. |4ll. This is consistent with the idea of extinction caused 
by absorption of dust associated with the GHIIR. Other dust-gas 
configurations would have caused a different set of maps. For ex- 
ample, if dust were behind the region, there would not have been 
a counterpart in the E{B — V) map to the peaks in the map at 
24 /im. In the right panel of Fig. |4] we present the absorbed Ha 
luminosity - defined as the difference between the total extinction- 
corrected Ha luminosity and the Hq luminosity co rrected for the 
foreground Galactic extinction, E{B - V)=0.044, iSchlegel et all 
[l998) - obtained using our derived reddening map with the 8 /xm 
contours overlaid. A comparison of these two panels shows how the 
24 nm emission presents a more compact distribution towards the 
centre of the region and correlates better with the absorbed Ha lu- 
minosity map than the 8 /im emission. This was seen in other H II 
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Fi gure 4. Left : Redd ening map for NGC 588 obtai ned assuming an intrinsic Ha/H/3 = 2.86 and a correction of 1 A in absorption for H/3. The extinction law 
of lpluks et al] h994) and E{B ~ V) = Ay/^-l faeke & LebofskJ 19851) were utilized. The map has been convolved with a Gaussian filter with cf=\" to 
better trace the extinction structure. The 24 ^m emission was overplotted with contours. The intensity contours are at 2, 5, 10, 20, 40, 60, 80, 95 per cent of 
the maximum intensity within the region. A 1% contour level coiTesponds to a 3(T value. Right: Absorbed Ha luminosity of NGC 588 with 8 /im emission 
contours overplotted. The contour levels coirespond to the same percentage as for the 24 /xm case. Here, a 1% contour level corresponds to a Ict value. 



regions like NGC 6 04 jRelano & Kennicutl|[2009l) and NGC 595 
( lReIanoetai]|201(]|) . The 24 /xm emission has proved to be a good 
tracer of the recent star formation for a large range of Hq luminosi- 
ties, rangi ng from moderate H II regions to dus ty powerful star- 
bursts (e.g. lRieke et al Jl 20091 : Icalzetti et alfeoiol) . The spatial cor- 
relation between the Hq luminosity and 24 /im maps indicates the 
infrared band as a star-formation tracer also in low dust environ- 
ments, such as NGC 588. Due to the low spatial resolution of the 
70 and 160 Spitzer bands we were not able to compare the 
emission at these wavelengths with our reddening map. Herschel 
data, which will be available soon to the scientific community, will 
permit a proper comparison to be made in the future. 

3.4 WR stellar population 

lRelanoetalJ ( l2010h presented a novel and simple technique to de- 
tect WR stars in a swift way and compared its results to classical, 
more time-consuming techniques. Using the same set of observa- 
tions (i.e. the datacube), one can identify the WR candidates by 
simulating the action of narrow filters and creating continuum sub- 
tracted maps at the emission bump at 4 700 A (the blue bump) and at 
5 700 A (the red bump), characteristic of WR stellar emission, and 
localizing the peaks of emission, afterwards. Then, the candidates 
can be confirmed by extracting the spectra of the associated spax- 
els. Following this methodology, we confirmed previous cataloged 
WR stars in NGC 595 and discovered a new one further away from 



the main ionizing clusters. Here, we apply the same methodology 
to NGC 588. 

There are two known WR stars in NGC 588 that have been 
widely studied in the past. The first one, name d UIT-011 by 
iMassev et all ( Il996l) - MC 3 bv lPrissen et al.l fiooj) - was detected 
for the first time using narrow band imaging by Iconti & Masseyl 
( 1981 ) and was spectroscopically confirmed by Massey & Contj 
( Il983t) later on. It was classified as WNL with Mv = -7.9 mag. 
The second one, named UIT-008 is a transition OfAVN9 star. 
Both sta rs were modelled using multi-band photometry with the 
HST bv lUbeda&DrissenI j200^ . They derived effective temper- 
atures of 57 000 K and 32000 K and bolometric luminosities of 
log(L/L0) = 6.48 and 5.97, for UIT-011 and UIT-008 respec- 
tively. 

Fig.|5]presents the continuum substracted blue bump map af- 
ter convolving with a 1" -Gaussian. We identify two main peaks 
of emission whose positions agr ee well with those previously re- 
ported for UIT-01 1 and UIT-008 jPrissen et alj2008l) . No aditional 
WR in NGC 588 was found. The extracted spectra for the two stars 
is presented in Fig.|6]and show clearly both the blue and red bump. 

The cases of NGC 595 and NGC 588, in H II regions, as well 
as existing ones for starburst galaxies (e.g. the Anntenae, II Zw 70, 
Mrk996, NGC 5253, IC 10 , Bastian et al. 2006; Kehrig et al. 200j; 
I James et al.l l2009l :l Monreal-Ibero et al.i.2010i : .Lopez-Sanchez et al.l 
2010) illustrate the effectiveness of IFS in the finding and charac- 
terization of the WR population. At this stage, the possibility of 
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Figure 5. 4650—4750 A map derived from the PMAS data after convolv- 
ing with a Gaussian of a=\'!Q. Continuum was subtracted by averaging the 
spectral ranges of 4490— 4540 A and 4755—4805 A. A logarithmic stretch 
covering ~0.3 dex was used to better enhance the peaks of emission. Labels 
with the detected WR are included and white contours delineate the spaxels 
utilized to create the extracted spectra. Contours correspond to the contin- 
uum subtracted R b road-band direct image from NOAO Science Archive 
iMassev et alj|2007l) . The orientation is north up and east to the left. The 
main ionizing cluster, at coordinates RA(J2000): Ih32m45.7s, Dec.(J2000): 
+30d38m55.1s, marks the origin of our coordinates system. 



using this technique routinely should be taken into account. In par- 
ticular, it would suit perfectly in the case of WR finding in galaxies 
at larger distances and, more important, with large gradients in their 
velocity fields. Here, the traditional technique of search for candi- 
dates via imaging first, and spectroscopic confirmation afterwards, 
might well miss some of the WR populations since the blue/red 
bump might move outside the spectral range of the narrow filter. On 
the contrary, the methodology presented here can be easily modi- 
fied and implemented to define what can be called synthetic tunable 
filters that take into account the movements of the galaxy and thus 
preventing these losses. An additional advantage of using an Inte- 
gral Field Unit, especially if it has a large field of view is the detec- 
tion of runaway WR sta rs ejected by the central star cluster of the 
region toravetalj200i) . In particular, IFS-based instruments with 
relativel y large field of view like PPak/PMAS ( Kelz et al. 2006), or 
MUSE (Baco n et aLlllOlOl) are (or will be soon) under operation. 
They open the possibility of carrying out surveys of large samples 
of galaxies where this kind of simple techniques could be particu- 
larly useful. 



3.5 Kinematics of tlie ionized gas 

The superior spectral resolution of the present observations, twice 
that of the observations of NGC 595, allowed us to derive the veloc- 



30 

g 25 

& 20 
= 15 

Li_ 

10 









UIT-01 1 












i - 


4000 


4500 


5000 5500 


6000 






\ 


UIT-008 


1 \ 











5000 

Uk 



Figure 6. Extracted spectra for two detected peak of emission in the PMAS 
continuum subtracted 4650—4750 A. Vertical arrows mark the position of 
the expected blue and red bumps. 
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Figure 7. Velocity field derived from the H/3 emission line. Contours cor- 
respond to th e continuum subtrac ted Ha direct image from NOAO Sci- 
ence Archive iMassev et al.ll2007l) . The orientation is north up and east to 
the left. The main ionizing cluster, at coordinates RA(J2000): Ih32m45.7s, 
Dec.(J2000): +30d38m55.1s, marks the origin of our coordinates system. 



ity field maps from the strongest emission lines. In Fig.|7]we show 
the map corresponding to H/3. No relevant differences were found 
from the map derived using [O III]A5007. The velocity field has a 
complex structure with values ranging between —190 km s~^ and 
— 110 km s^^. The north-west part of the region seems to be more 
redshifted than the south-east, which is the region with higher H/3 
surface brightness (see Fig.[3]l. In the surroundings of the location 
of the stellar cluster, marked in Fig. |7] as a black cross, there is a 
pronounced velocity gradient: the north-west part has velocities of 



NGC 588 with PMAS 9 



— 120 km s^^ while the south-east has velocities of —170 km s~^. 
The velocity separation of --^25 km between these zones and 
the location of the stellar cluster, corresponding to a velocity of 
^-^30 km s ~^ in the galaxy plane (inclination of M 33, i = 56 deg 
Ivan den B ergh 2000), and the symmetry of the velocity field sug- 
gest the existence of a shell expanding in the interior of the region. 
The shell expansion velocity, ~30 km s~^, is slightly lower than 
the values observed in high luminos ity H I I regions of a set of spiral 
galaxies (Vexp ~ 40 — 90 k m s~^. | Relano et al . 2005) and also in 
NGC604(Ve^p -^40 km s~ MYang et al.ll993). However, based on 
a kine matic study of NGC 588 and NGC 604. lMunoz-Tunon et al] 
( Il996h suggest that NGC 588 is more evolved than NGC 604 and 
therefore, we would expect lower velocities for the shells in the first 
region than in NGC 604, consistent with the result found here. 

In order to check whether the winds coming from the stellar 
population within the H II region could produce the expansion of 
the observed shell, we have made a crude estimation of the kinetic 
energy involved in the shell and compared with the input kinetic 
energy from the stars. Usin g the Ha luminosity of the region from 
iRelano & Kennicut^ j2009l) , we predict an emission measure (EM) 
of 4000 (pc cm~®) for an H II region radius of 140 pc, correspond- 
ing to the aperture radius used to obtained the Ha luminosity. The 
EM is then used to derive a < rie >rms of 5 cm"'' and integrat- 
ing over the H II region volume we derive a total ionized mass of 
~6x 10''' Mq. Assuming, as an upper limit that the whole mass has 
been swept up by the shell we obtain a kinetic energy for the shell 
of 3.8 X 10^^ erg. Starburst99 models (Leitherer et al. 1999) using 
ranges of values for the stellar mass of 1 - 6 X 10^ Mq and age 
of 3.5 — 4.2 Myr, a Salpeter initial mass function, and metallici- 
ties of Z=0.004 and 0.008 give a range of kinetic energy input of 
5.9 — 64.2 X 10^^ erg, at least twice as high as the upper limit of 
the kinetic energy of the shell. This crude calculation shows that 
the winds from the stellar cluster within NGC 588 are able to pro- 
duce the observed shell in the H II region and create the observed 
H/3 morphology with holes and filaments shown in Fig. (3) More- 
over, the lower energy associated with the shell, 30 km s~^, com- 
pared w ith values closer to 100 km s"^ for younger, less evolved 
regions ([Relario et al.ll2005l) . is consistent with with the findings by 
iMufioz-TuSon et al.Ul99^ . 



3.6 Characterization of the ionized gas 

3. 6. 1 Density structure 

Electron density can be determined from the ratio between two 
lines of the same ion emitted by different levels with similar ex- 
citation energies. We used the [S II]A6717/[S II]A6731 ratio in our 
analysis. As is shown in the lower right map of Fig. [8] no structure 
for the Tie was found. We measured a mean(± standard deviation) 
[S II]A6717/[S II] A6731 value of 1.20(±0.17), which is consistent 
with the value derived from the integrated spectrum (see Table [2l(. 
For the assumed temperature, this implies a Ue of ~250 cm""^ 
and agrees w i thin t he un certainties with th e values reported by 
IVflchez et all h988l) and IJamet et alj ( l2005h . This value for the 
electron density corresponds to the density of the clumps within 
the region and differs from the r.m.s. electron density derived in the 
previous section using the Ha surface brightness of the region. The 
ratio of both density estimates is a measure of the volume fraction 
occupied by dense clumps. 



3.6.2 Line ratios in the BPT diagnostic diagrams 

Diagnostic diagrams, where different areas of a given diagram 
are occupied by gas excited via different mechanisms, have been 
widely used to study the ionization conditions of the ISM. In the 
optical spectral range, th e most popular are pr obably the BPT dia- 
gra ms, first proposed by Baldwin et al.l ( Il98ll) and later reviewed 
by IVeilleux & Osterbrockl jl987l) . Their wide use to study the 
ionization conditions in sta r-forming and starburst galaxies (e.g. 
lAlonso-Herrero et al.l l201(ih is due to the fact that they involve 
emission lines that are relatively strong and line ratios that have 
(almost) no dependence on the extinction. 

With the present 2D unbiased mapping, one can determine 
these line ratios locally. In this way, it is possible to evaluate their 
dependence on the position within the GHIIR and relative surface 
brightness of the area under study. Moreover, one can make a com- 
parison between integrated and local values. 

We present the maps for the three available line ra- 
tios involved in the BPT diagrams - namely [N II] A6584/Ha, 
[Sll]AA6717,6731/Ha, and [O IIl]A5007/H/3 - in Fig. [8] These 
maps show that the ionization structure in NGC 588 is complex. 
The [Nll]A6584/Ha and [S II] AA6717,6731/Ha maps present a 
rather similar structure. In both cases, the minimum is located nei- 
ther at the peak of emission in H/3 (i.e. ionized gas) nor at the one 
for the continuum (i.e. stars) but in the middle point between them. 
Then, line ratios increase outwards, following the ring structure of 
the region. 

The [OIII]A5007/H^ map is roughly similar to the 
[Nll]A6584/Ha and [S II] AA6717,6731/Ha maps but its value 
varies in the opposite sense. Moreover, there are two main 
differences. The peak of [O III] A5007/H/3 (maximum for this 
line ratio) is broader than the peak for [N II] A6584/Ha and 
[S II]AA6717,6731/Ha (minimum for these ones). More relevant, 
there is an area of ~ 5" x 13" centred at ~ [7'.'0, -S'.'O] of 
elevated [O IIl]A5007/H/3 values that do not show specially low 
[Nll]A6584/Haand [S II]AA6717,6731/Ha. 

In order to better assess how the ionization conditions change 
in different parts of NGC 588, we divided our data in three lu- 
minosity bins, which sample the low, medium and high surface 
brightness areas of this GHIIR. In Fig. |9l we present the posi- 
tion of each individual spaxel in the BPT diagnostic diagrams to- 
gether with the borders that separate H II region-like ionization 
from ionization by other mechanisms according to several authors 
jVeilleux & Osterbrock'l987l: lKewlev et alj200lljKauffmann et al.l 
12003 : Stasiriska et al. 2006). As expected, all line ratios are within 
the typical values expected for an H II region-like ionization. 

However, there are differences between the diagrams associ- 
ated with the different luminosity bins. Firstly, the range of ob- 
served values varies from ~0.3 to ~0.6 dex for [O IIl]A5007/H/3, 
from ~0.6 to ~I.I dex for [N II]A6584/Ha and from ~0.5 to 
~I.I dex for [S II]AA6717,6731/Ha, with larger ranges in those 
areas with lower surface brightness Secondly, as indicated by 
the position of the red triangles, higher [N II] A6584/Ha and 
[S II] AA6717,6731/Ha ratios and lower [O IIl]A5007/H/3 ratios are 
detected in the areas of lower surface brightness. This implies that 
the degree of ionization gets smaller with increasing distance from 
the ionizing source. 

How do these results compare with our findings for NGC 595? 
In general, the tendencies in the differences between the integrated 
values and those for the individual spaxels in NGC 588 are similar 
to those found in NGC 595. However, the range of observed val- 
ues in NGC 588 is smaller. This is seen in [S II]AA6717,6731/Ha 
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Figure 8. Map of the observed [N II] A6584/Ha (upper left), [S II]AA6717,6731/Ho (upper right), [O III]A5007/H/3 (bottom left) and [S II]A6717/[S II] A6731 
(bottom right) line ratios. Contours coirespond to the continuum subtracted Ha direct image from NOAO Science Archive (,Massev et al. 2007). The orientation 
is north up and east to the left. The main ionizing cluster, at coordinates RA(J2000): Ih32m45.7s, Dec.(J2000): +30d38m55.1s, marks the origin of our 
coordinate system. 
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Figure 9. Position of the individual spaxels in NGC 588 in the BPT diagnostic diagrams separated in thre e h(R0)/L(Hl3)max bi n s (Left column: 0.00-0.25; 
Middle column: 0.25-0.50; Right column: 0.50-1.00). Solid curves show the empirical borde rs found bv |yeilleux & Osterbrockl j 19871) between ionization 
caused by different mechanisms, while dotted lines show the theoretical borders proposed bv lKewlev et alj )200ll) to delimit t he ar ea where the line ratio s 
can be explained by star formation. Black dashed and dot-dashed lines show the revised borders bv iKauffmann et al] j2003|) and IStasiilska et si\ j20O6l) . 
respectively. These were empi rically determined u sing SLOAN data. Green diamonds mark the values reported bv lVflchez et"d] <1988l) ^ while yellow inverted 
triangles are used for those of 'jamet et ajj j2005t) . Red triangles and blue squares indicate the mean values for a given bin and the values derived from the 
integrated spectrum as reported in Table|2] respectively. The number of considered data points as well as the L(H/3)/L(H/3)maa; range are indicated in the right 
upper comer of the diagrams involving [N II] A6584/Hq. Typical errors are shown in the left lower comer of each diagram. 
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Figure 10. Maps of the hne ratios sensitive to the ionization parameter: [O III] AA4959.5007/[ OlI] AA3726,3729 (le ft) and [O II] AA3726,3729/H/3 (right). 
Contours con'espond to the continuum subtracted Ha direct image from NOAO Science Archive I'M assev et al.ll2007l) . The orientation is north up and east to 
the left. The main ionizing cluster, at coordinates RA(J2000); Ih32m45.7s, Dec.(J2000): -H30d38m55.1s, marks the origin of our coordinates system. 
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Figure 11. Diagrams involving two U sensitive line ratios. Upper row: [O III]AA4959,5007/[OlI]AA3726,3729 vs. [S II] AA6717,6731/Hq. Lower row: 
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and specially in [O IIl]A5007/H/3, where the range for NGC 595 
is twice as large and can be understood in terms of the different 
mapped area. While in NGC 588, we map just up to the border of 
the ring, in NGC 595, we were able to go into the very low surface 
brightness component, further away from the shell. 

An interesting result stands out in Fig. |9] after comparing the 
mean line ratios (red triangles) for each bin with those measured 
for the integrated spectrum (blue squares): the case of NGC 588 
shows that independently of the utilized line ratio, integrated val- 
ues are more representative of the ionization conditions in the low 
surface brightness areas (L(H/3)< 25%L(H/3)„iai), which occupy 
~70% of the region, than of those in the brightest parts. Moreover, 
[Nll]A6584/Ha and [S II]AA6717,6731/Hq ratios derived for the 
integrated spectrum differ by ~0.3 dex from those derived with 
long-slit (i.e. green diamonds and yellow inverted triangles). Sim- 
ilar effects have been fo und in the few GHIIRs mapped up-to-dat e 
at such a level of detail jPellegrini et"ai]|2010l : iRelano et alj|2010h . 
However, when observing H II regions with long-slits, these are 
usually the selected areas to be observed. Thus, this result should 
be taken into account when interpreting the ionization conditions 
in GHIIRs in distant star-forming galaxies, specially if calibrations 
derived from observations of local H II regions are utilized. As an 
example, at ~ 25 — 40 Mpc, NGC 588 would occupy ~1 arcsec^ 
on sky, which is the typical size that nowadays IFS-based instru- 
ments can resolve under typical seeing conditions. At further dis- 
tances, the situation would become even more uncertain since a 
typical spaxel would sample in addition some emission associated 
with the Diffuse Ionized Gas. 

3.6.3 Ionization parameter 

Line ratios presented in previous section are useful to study the 
ionization properties of a given region/galaxy. However, to have 
a more detailed view of the ionization structure, one should 



estimate how physical-chemical quantities like metallicity, relative 
abundances and ionization parameter vary within the region. In 
this section, we will explore the behaviour of the different tracers 
of the ionization parameter {U). This can be estimated from 
the ratio of lines of the same element that trace two different 
ionization states (e.g. [O IIl]AA4959,5007/[O II] AA3726,3729). 
Assuming that the metallicity is known, one can also use 
the [OlIl]A5007/H/3 or the [S II]AA6717,6731/Hq line ra- 
tios. The [Sll]AA6717,6731/Ha is presented in Fig. [H while 
the maps for [O III] AA4959,5007/[O II] AA3726,3729 and 
[Oll]AA3726,3729/H/3 appear in Fig. [TO] In all three cases the 
observed structure is the same: ratio values corresponding to high 
ionization parameters are found between the peak of emission in 
H/? and the main ionizing cluster, while ratio values typical of 
lower ionization parameters are found outwards, following the ring 
structure of the nebula. 

This is better seen in Fig. [TT] which shows a good correla- 
tion between the different tracers for the three flux bins under 
consideration. Also, as in the BPT diagrams, integrated line ra- 
tios are not dominated by the brightest zones of the gas but by 
the larger low surface brightness areas. In particular, differences in 
the [O III] AA4959,5007/[O II] AA3726,3729 ratio between the inte- 
grated values and those of the spaxels with high surface brightness 
can be of ~ 0.5 dex on average and as high as ^ 0.8 dex. 

Do these line ratios make consistent predictions of the 
ionization parameter? For the purpose of this di sc ussion we 
will use the expressions provided by iDiaz et ai] ( I2OOOI) as- 
suming a metallicity of 0.3 Zglfl Note that small variations 
of the metallicity (i.e. allowing for a range between 0.2 and 
0.4 Zq) would imply an offset in the estimated logL'^ be- 



3 We have employed 12 + log(O/H)0 = 8.66, from lAsplund et al.l j2004h 
and the metallicity derived in Sec. 13. II 
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tween -0.20 and 0.15 dex. Also, in order to minimize the ef- 
fect of the extinction, we utilized the [S II]AA6717,6731/Ha 
line ratio instead of the [S Il]AA6717,6731/H/3 ratio and as- 
sumed H/3=Ha/2.86. Fig. [TT] also contains the locus of line ra- 
tios that trace the same ionization parameter and show that accord- 
ing to these relations, [O IIl]AA4959,5007/[O II] AA3726,3729 and 
[O II] AA3726,3729/H/3 predict relatively consistent results while 
the [S II] AA67 17,673 1/Ha would correspond to smaller ionization 
parameters, even when considering the integrated spectrum. In any 
case, in the quantification of the ionization parameter we are not 
taken into account the fraction of i onizing phot ons leakin g the H II 
regio n which can be up to ~50% jZurita et~ar,,2000 ; Relano et al.l 
l2002h . 

Both, the spatial variations of the different line ratios (see 
Fig.[8]l and the observed excess in the [S II] AA67 17,673 1/Ha ratio 
when compared with photoionization models (see Fig. II It are a di- 
rect consequence of the ionization structure of the GHIIR and con- 
stitute a nice observational counterpart to the 3D modelled struc- 
ture of ionized regions. These changes are obvious across the maps 
whereas the comparison of the different line ratios in Fig.ll 1 [traces 
the ionization structure in the line of sight since areas of differ- 
ent degree of ionization ar e traced by differe nt ions. Specifically, 
the expressions provided bv lDfazetal.] ( l2000l) were derived for the 
integrated spectra produced by one single ionizing star. However, 
the situation in a GHIIR like NGC 588, where the ionizing stars 
are distributed in 3D in an irregular manner, is much more com- 
plex. In this scheme, the characteristic size of the different zones 
of the ionization structure will be determined by the architecture of 
the GHIIR (i.e. by the relative distribution of the ionizing sources). 
Thus the lower ionization species such as 5^, will delineate the 
more extended and common component while O'^^ will be con- 
fined to different high ionization zones at the vicinity of the ion- 
izing sources. A model predicting the two-dimensional observable 
structure of the region designed to match our observations is in 
preparation (Perez-Montero et al. in prep.). 

Using the [O IIl]AA4959,5007/[O II]AA3726,3729 as a base- 
line, we can comp are our results with those for NGC 595 
jRelano et al] |2010|) . NGC 588 presents higher values of 
[OIII]AA4959,5007/[OII]AA3726,3729, which would imply dif- 
ferences in U ranging between -^0.7 and ^^0.4 dex, being these dif- 
ferences larger when we are closer to the main ionizing cluster. At 
similar gas densities (as is the case for NGC 595 and NGC 588), the 
ionization parameter depends on the characteristics of the ionizing 
stars, the geometry of the region and the filling factor. Regarding 
the stars, the hotter these are, the higher number of ionizing pho- 
tons they produce, and thus, a higher ionization parameter is locally 
expected. In general, the lower the metallicity of the region is and 
the younger the stellar population is, the larger number of hot stars 
is expected. These tendencies can be seen by modelling of inte- 
grated spectra of H II regions ionized by given stellar populations 
(e.g.^Levesque et al. 2010). Thus, the lower metallicity of NGC 588 
and the youth of its stellar population with respect to NGC 595 
can, at least partially, explain the difference between the observed 
[OIII]AA4959,5007/[OII]AA3726,3729 line ratios. However, 2D 
detailed modelling also showed that the relative distribution of the 
ionizing sources is an i mportant parameter (Ercolano et al. 2007; 
Ijamet & Morissetll2008l) : GHIIRs with more sparsely distributed 
ionizing sources have lower ionization parameter. In that sense, a 
detailed modelling of the region will help to disentangle the relative 
role of geometry and filling factor (Perez-Montero et al. in prep.). 



3.6.4 Metallicity tracers 

Ideally, metallicity is calculated in a direct manner. This requires 
the determination of the electron temperature via detection of the 
- e.g. - faint [OlIl]A4363 line. Another possibility is the use of 
certain combinations of strong emission lines for which empirical 
and/or theoretical calibrations have been established. Here, we will 
focus on those that can be evaluated using emission lines within our 
spectral range. In particular, we will see the spatial distribution of 
the metallicity tracers as well as explore their reliability as proxies 
of the metallicity. 

The map for N2 = log([N II]A6584/Hq) was presented in 
the upper left corner of Fig. [8] while those for the J?23 and the 
N203 parameters appear in Fig.[T2] None of them presents a uni- 
form distribution. A comparison of the maps presented in Figs. 
[T2I and [To] shows how the N203 and 7V2 parameters are mod- 
ulated by U (i.e. the higher values of U correspond to the ar- 
eas showing the lower A'^203 and N2 values). The R23 parame- 
ter presents a different distribution due, in great measure, to the 
anomalous high [O III] A5007/H/3 values measured at the south- 
east of NGC 588 (see section |3^.6.2| l. This quadrant present typical 
values of log(i?23)^0.9 very much at the turnover region of the 
i?23 — Z relation while the other areas of NGC 588 have a rela- 
tively uniform distribution with typical values of log(i?23) ~0.8- 
0.9. As we stated in Sec. 13.11 given the measured values and that 
the Z — R23 relation is two-valued, R23 is not a reliable metal- 
licity tracer for NGC 588 and will not be discussed further in this 
section. 

The dependence of the remaining metallicity tracers on 
the ionization parameter is better seen in Fig. [13] where we 
grouped the data in three bins of degree of ionization using 
[O III] AA4959,5007/[O II]AA3726,3729 tracer as basefine. iV203 
varies from —2.8 to —1.3 dex when going from high to low val- 
ues of the ionization parameter while A^2 varies from —1.9 to 
—0.8 dex. Also, the comparison of the mean and integrated ratios 
shows how measurements of N203 and N2 are dominated by the 
areas with low U, which mostly correspond to zones of low surface 
brightness. 

How these variations in the measured ratios translate 
into uncertainties in the metallicity determination? To an- 
swer this question, we used the calibrations proposed by 
jPerez-Montero & Continil j2009[) for N2 and iV203. The relation 
between the predictions of the two parameters is shown in Fig. 1141 
The integrated spectrum gives consistent estimates of the metallic- 
ity for both tracers with a mean of 12+\og{0 / H} — 8.16 (see Sec. 
|3.1b . and in agreement with those reported bv ljamet et alj ( l2005h . 
Since the utilized calibrations were derived empirically using data 
for integrated spectra, this is not particularly unexpected. Hereafter, 
this value will represent the reference metallicity for NGC 588. 

The first result extracted from Fig. [14] is that independently 
of the utilized tracer, we derive different values for the metallic- 
ity which depend on the position in the H II region at which the 
measurement is taken: the range of predicted metallicities covers 
~0.6 dex when using N203 and up to ~1.0 dex for N2. Thus, 
N203 seems to be more reliable metallicity tracer than N2. 

The second result is that the values of the metallicities pre- 
dicted in the zones of low excitation using the different in- 
dices show less variation between them than those in the high 
excitation zones. Both indices predict rather low metallicities 
in the latter, but the effect is bigger for A'^2, where the rel- 
ative underestimate can reach ~0.5 dex. The ionization struc- 
ture of the region gives a good explanation for this. Where the 
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Figure 12. Maps of the different observed line ratios that trace the met allicity: left: th e R23 parameter; right: the N203 parameter Contours con'espond to the 
continuum subtracted Hq direct image from NOAO Science Archive i Massev et al . 2007). The orientation is north up and east to the left. The main ionizing 
cluster, at coordinates RA(J2000): Ih32m45.7s, Dec.(J2000): +30d38m55.1s, marks the origin of our coordinates system. 
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Figure 13. [Nii]A6584/Hq! vs. 7V203. The color/symbol code is as in Fig. |9] The number of data points as well as the range of 
log([0 III] AA4959,5007/[O II]AA3726,3729) considered are indicated in the lower right comer of the individual diagrams. 



[OIII]AA4959,5007/[OII]AA3726,3729 ratios are low only gas 
with a relatively low degree of ionization emits along the line of 
sight, but where higher values of this ratio are found there is a sig- 
nificant contribution from zones with a higher degree of ionization, 
closer to the ionizing stars. This leads to a d eficiency in the [N II] 
line intensity (see e.g. iLevesque et"ai]|201(]|) which leads to an un- 
derestimate of the metallicity when N2 is used as the tracer. 



3.6.5 Relative abundance tracers 

Massive stars, via their strong winds and supemovae explo- 
sions, eject all the heavy elements previously synthesized con- 
tributing to the metal-enrichment of the gas. In particular, dur- 
ing their main sequence to earlier WR phases, material pro- 



cessed through the CNO cycle and enriched with nitrogen is 
poured into the ISM. This contamination can be observationally 
spotted by the detection of areas in the ISM with higher rela- 
tive abundance in nitrogen. The chances for these detections are 
a priori low, since the time scales involved in the dilution of 
this nitrogen are relatively quick (e.g. Monreal-Ibero et al. 201(i)- 
However, this kind of data, that maps the ISM from the clos- 
est vicinity of the polluting stars to the more external parts of 
the GHIIR offer an invaluable opportunity to look for such in- 
homogeneities. Here we will look for areas of enhanced nitrogen 
abundance by means of the [N II]A6584/[S II] AA6717,6731 and 
[Nll]A6584/[0 II] AA3726,3729 ratios. 

The maps for our two relative abundance tracers under 
consideration are displayed in Fig. [15] No inhomogeneities 
are detected in any of them, not even in the vicinity of 
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Figure 14. Compari son of the metallicity e stimates from N203 and N2. The color/symbol code is as in Fig. [9] however we include the metalUcities from 
IVilchez et al]h98g|) and ljamet et aiWim f) derived using the direct metliod instead of those from the and N203 line ratios. The red bars associated with 
the me an values indicate the standard deviation. In addition, the expected metallicity from the metallicity gradient for M 33 according to Rosolowsky & Simq^ 
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in the lower right comer of the individual diagrams. The locus of equal estimated metallicities is indicated with a black line. The data point distributions as 
well as their fit to a Gaussian are shown in the margins. 
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Figure 15. Maps for line ratios tracing inhomogeneities in the abundances. Left: [N II] A6584/[ S II] AA67 17,673 1. R ight: [N II] A6584/[0 II] AA3726,3729. 
Contours correspond to the continuum subtracted Ha direct image from NOAO Science Archive i Massev et alj|2007h . The orientation is north up and east to 
the left. The main ionizing cluster, at coordinates RA(J2000): Ih32m45.7s, Dec.(J2000): +30d38m55.1s, marks the origin of our coordinates system. 



the two WR stars. The distributions for both line ratios (not 
shown) can be well reproduced by a single Gaussian func- 
tion with mean(±standard deviation) of — 0.34(±0.08) and 
-1.27(±0.10) for the log([N II] A6584/[S ii]AA6717,6731) and 
los([Nll]A6584/[ Oll]AA3726.3729). respec t ively. Using the cal- 
ibrations given by IPerez-lVIontero & Continil ^20091) . these values 



imply a \og{N/0) — —1.29 ± 0.10. This is in good agreement, 
within the uncertainties with the nitrogen abundance reported by 
Ijamet et al.j ^2005) and higher by a factor of ~1.7 to the one re- 
Dorted bv lVilchez et al.l ( Il988l) . 
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4 SUMMARY AND CONCLUSIONS 

This work presents a detailed analysis of NGC 588, a GHIIR 
in M33. The study is based on the joint analysis of more than 
1 000 optical spectra mapping in an un-biased manner an area of 
~120 pcxl80 pc, together with mid-IR images in the 8 /im and 
24 /xm Spitzer bands. This allowed us to better understand the re- 
lationship between the different elements (i.e. stellar populations, 
ionized gas and dust) playing a role in the region. Our main results 
are: 

1 . Flux averaged properties of the region were derived by us- 
ing the integrated spectrum. Differences in the relative line inten- 
sities from those previously reported (and usually sampling the 
brightest area of the region) range between the 8% and the 50%. 
The He derived from the [S II] line ratio is consistent with being 
within the low density limit, as reported in previous works. Our 
derived metallicity is in agreement with previous measurements as 
well as with the expected metallicity according the metallicity gra- 
dient for M 33. 

2. This two-dimensional method reveals complex structure in 
the extinction distribution which had not been found with previ- 
ous long-slit measurements. Our optical reddening map presents 
three maxima that correlate well with those for the dust emission 
at 24 /im and 8 /im. Moreover, the absorbed Ha luminosity map 
reproduces the structure observed in the 24 fim image from Spitzer 
supporting this band as a good tracer of recent star formation. 

3. Using the same methodology as in iRelanoetalltold) , we 
confirm the location of the two already detected WR stars in this 
region. This consolidates this technique as an efficient manner to 
find WRs. However, no new WR was detected. Since we mapped 
almost its whole surface, the census of the WR stellar content of 
NGC 588 can be considered complete. 

4. We derived velocity maps from the strongest emission lines. 
No remarkable differences were found between them. We measured 
a velocity difference of ~ 50 km s~^ between the areas of high and 
low surface brightness. This implies a smaller kinetic energy for the 
expanding shells than the one expected from the stellar population 
of NGC 588. The velocity field is consistent with NGC 588 being 
an evolved H II region, in agreement with previous studies. 

5. We measured the electron density over the whole face of the 
region using the [S II] emission line ratio and no significant vari- 
ations for Tie were found. The estimated values for the [S II] ratio 
are within the low density regime with a mean value corresponding 
to 7ie=250 cm~'^ for the assumed temperature. This is in agreement 
within the uncertainties with the values reported for specific areas. 
Using the EM of the region reported in the literature, we estimate a 
< ne > of 5 cm~^ for the whole GHIIR. 

6. [Nii]A6584/Hq and [S II] AA67I7,673I/Ha present a ra- 
dial structure with values increasing from the centre towards the 
outer parts of the region. In general, [O III] A5007/H/3 shows an 
inverted pattern. However, we detected an area towards the south- 
east of the region with high [O III] A5007/H/3 line ratios but not 
particularly low [N II]A6584/Ha and [S II]AA6717,6731/Ha ra- 
tios. Line ratios involved in the BPT diagrams are typical of pho- 
toionization caused by massive stars. However, a larger range of 
observed line ratios is found at lower surface brightness, vary- 
ing from ~0.5 to ~1.2 dex for [N II]A6584/Ha, from ~0.7 to 
~1.7 dex for [S II] AA6717,6731/Hq, and from ~0.3 to ~0.5 dex 
for [OIII]A5007/H^. 

7. We studied the behaviour of three tracers of the ioniza- 
tion parameter. In all three cases the ratios corresponding to high 
ionization parameter are found between the peak of the emis- 



sion in H/3 and the main ionizing source decreasing when go- 
ing outwards. Also, differences between the integrated and lo- 
cal values of the U tracers can be as high as ~0.8 dex, spe- 
cially when using the [O IIl]AA4959,5007/[O II]AA3726,3729 ra- 
tio and in the high surface brightness spaxels. [O III] A5007/H/3 and 
[OIII]AA4959,5007/[OII]AA3726,3729 predict similar local ion- 
ization parameters and consistent with those estimated for the inte- 
grated spectrum of an H II region ionized by a single star. However, 
the [S II] AA67I7,673I/Ha line ratio departs from these predictions 
reflecting the more complex ionization structure in GHIIRs. The 
higher [O III] AA4959,5007/[O II] AA3726,3729 ratios in NGC 588 
than those found in NGC 595, specially close to the main ioniza- 
tion source, can be partially explained by the lowest metallicity of 
NGC 588 and the youth of its stellar population. 

8. -R23 presents a relatively uniform distribution, with varia- 
tions within ~ 0.2 dex. However, since it lies at the turnover of 
the R23-Z relation, R23 cannot be considered a good metallicity 
tracer for this particular region. Other metallicity tracers are mod- 
ulated by the degree of ionization. The range of measured values 
for both the N203 and the tracers is ~1.5 dex. Both ratios 
vary with U with variations of up to ~0.5 dex on average when 
going from high to low ionization parameter. This implies an un- 
certainty in the determination of the metallicity associated with the 
position where the measurement is made of ~0.6 dex and ~I.O dex 
for N203 and N2, respectively. Therefore, N203 seems to be a 
more accurate metallicity tracer than A'^2. 

9. The N/O relative abundance is homo- 
geneous over the whole face of the region ac- 
cording to the [Nll]A6584/[Sll]AA6717,6731 and 
[Nll]A6584/[Oll]AA3726,3729 line ratios with typical val- 
ues of log(iV/0) = -1.29 ±0.10. 

10. A common result emerges from our analysis of the line 
ratios involved in the BPT diagrams and those tracing metallicity 
and ionization parameter: the comparison between local and in- 
tegrated values shows that the line ratios from GHIIR in galax- 
ies at distances ^25 Mpc are likely to be dominated by the 
ionization conditions in their low surface brightness areas (i.e. 
L(H/3) < 25 %h(lil3}m ax) ■ This agrees with the results found by 
|Pellegrinietai] l l20IOl) after analyzing 2D line ratio maps of 30 Do- 
radus at the LMC. 
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